Photodynamics of stress in clamped nematic elastomers 
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We describe the complex time-dependence of the build up of force exerted by a clamped photo- 
elastomer under illumination. Nonlinear (non-Beer) absorption leads to a bleaching wave of a 
significant cis isomer dye concentration deeply penetrating the solid with a highly characteristic 
dynamics. We fit our experimental response at one temperature to obtain material parameters. 
Force-time data can be matched at all other temperatures with no fitting required - our model pro- 
vides a universal description of this unusual dynamics. The description is unambiguous since these 
are clamped systems where gross polymer motion is suppressed as a possible source of anomalous 
dynamics. Future experiments are suggested. 

PACS numbers: 83.30.Va, 61.41. +e, 61.30.-v, 78.20.H-, 81.40.Jj, 82.50.Hp 



I. INTRODUCTION 

Unusual properties of liquid crystal elastomers (LCE) 
arise from a coupling between the liquid crystalline or- 
dering of mesogenic molecules and the elasticity of the 
underlying polymer network. Crosslinked networks of 
polymer chains of a LCE include mesogenic units which 
either belong to the polymer backbone (main-chain LCE) 
or to side units pendent to the backbone (side-chain 
LCE) 1]. The shape of a monodomain nematic LCE 
strongly depends on the temperature dependent nematic 
order parameter Q(T). This connection is a consequence 
of the coupling of Q(T) with the average polymer chain 
anisotropy. By increasing the temperature of a LCE, 
Q(T) decreases which leads to a decrease of the polymer 
backbone anisotropy. This decrease of anisotropy mani- 
fests itself as an uniaxial contraction of the LCE sample. 

Mechanical change can be realized also through a 
change of nematic order by other means. As first shown 
in LCEs by Finkelmann et al. @, changes can also be 
achieved by introducing photo-isomerizable groups (e.g. 
azobenzene) into their chemical structure. These struc- 
tures will be referred to as nematic photo-elastomers. By 
absorbing a photon, azobenzene dye molecules can make 
transitions, with a quantum efficiency rj t , from their lin- 
ear (trans) ground state to the excited bent-shaped (cis) 
state. While the rodlike trans molecules contribute to 
the overall nematic order, the bent cis molecules act as 
impurities which reduce the nematic order parameter and 
lower the nematic-isotropic transition temperature. The 
illumination of photo-elastomers causes the reduction of 
nematic order, which in turn produces an uniaxial con- 
traction of the sample. On switching off the irradiation, 
the nematic order parameter recovers its initial dark state 
value, which results in a macroscopic expansion of the 
sample. 
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The characteristic time of the mechanical response 
of the sample is in the rather wide range of millisec- 
onds [1, H[ to hours @, To see whether slow on- 
response in Q was caused by a slow polymer dynamics of 
the nematic elastomer, or by a slow photo-isomerization 
kinetics, experiments were done [5- 7] on clamped photo- 
elastomers. Contrary to the case of the measurements on 
a freely suspended sample, the clamped setting requires 
no physical movement of chains in the network. In such 
settings the sample is fixed at both ends of its longest 
dimension (along the nematic director). Then, upon ir- 
radiation the sample cannot shrink and a retractive force 
is exerted on the clamps - see Fig. Q] for an illustration. 
This force is measured as a function of time until the 
photo-stationary state is reached. On switching off the 
irradiation the relaxation of the force begins. It was first 
observed by Cviklinski et al. @ that the dynamics of the 
nematic order parameter matches the dynamics of the 
mechanical response meaning that the rate-limiting pro- 
cess is dominantly the photo-isomerization. They also 
found that their systems display simple exponential pro- 
cesses for the build up of and decay of force. 

Subsequently, it has been shown in other systems Q 
that the stress - temperature experimental data for 
the on-process can be approximately fitted to a simple 
stretched exponential form a n (t) = 1 — exp[— (t/r n) ]■> 
where a n (t) is the normalized stress exerted on the 
clamps (stress at time t divided by the photostation- 
ary stress), and r on and /3 on < 1 are fit parameters. 
Similarly, it was found that the normalized stress in 
the off-process can be fitted to the nearly exponential 
law a n (t) = exp[— (t/r ff) /3off ], where f3 a g w 0.9. Such 
complex dynamical response cannot be attributed (as is 
usual) to polymer dynamics since we have clamping. In 
this paper we show that these findings can in fact be 
successfully described by our simple model of photody- 
namics and its conversion into stress in clamped nematic 
elastomers. 

According to the Beer law of light absorption the light 
propagating in a thick absorbing sample is attenuated at 
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constant rate. The light intensity at a depth x into the 
sample is I{x) — I e~ x / dt , where Iq is the incident inten- 
sity, and dt is the characteristic penetration depth of a 
given material due to absorption by trans isomers. How- 
ever, it has been shown that the simple Beer law for light 
attenuation through the sample containing dye molecules 
might be inaccurate, due to the so-called photobleaching 
effect [SrllH. This effect is caused by depletion of trans 
isomers, which allows light to penetrate to greater depths 
than those predicted by Beer's law. If dye molecules that 
absorbed photons don't return to their trans state im- 
mediately, the new photons falling on the sample cannot 
be absorbed in the initial layers and therefore propagate 
through the sample following a nonlinear absorption law. 

Using the model described in Sec. HH we calculate the 
stress exerted on the clamps during and after light irradi- 
ation. Then we relate the stress to the light absorbancc 
A = — ln(J/7o) at the back of the sample. We apply 
the full nonlinear absorption model which takes into ac- 
count the forward trans to cis and back cis to trans 
photo-isomerization as well as thermal-isomerization of 
cis molecules back to the trans state [llj ]. The major 
puzzle in photo-actuation is thereby addressed. Beer 
penetration depths are typically in the 1 — 10 /im range 
for normal dye loadings and hence orders of magnitude 
less than the sample thickness. If a small volume frac- 
tion of solid is photo-contracted, one expects the overall 
mechanical response to be small in the Beer limit. We 
show that the stress is proportional to the cis concen- 
tration and thus bleaching allows an appreciable sample 
volume to contribute to the force as a wave of trans to 
cis conversion deeply penetrates. To test the validity of 
our model we fit the on-process a n (t) data of Ref. |8| at 
one temperature in order to fix the material constants 
determining the photo-processes. The W n (t) response 
at other temperatures then does not need fitting - the 
same material constants suffice, after they are shifted by 
the separately measured changes in thermal relaxation 
times with temperature change, to reproduce the stress 
response. This remarkable universal agreement between 
theory and experiment confirms the hypothesis of the 
domination of photo-isomerization over polymer dynam- 
ics. We hence explain the observed stretched exponential 
(fion < 1) processes in terms of nonlinear spatio-temporal 
photodynamics, rather than the usual (unknown) poly- 
mer relaxation processes normally assumed to be behind 
such complex dynamics. Our analysis reveals that it is 
important to consider back cis to trans photoconversion. 



II. MODEL 

In this Section we consider a simple model of stress 
dynamics of clamped nematic photo-elastomers. The re- 
sulting stress is calculated within a nematic rubber model 
in Section III Al while the process of nonlinear light ab- 
sorption is outlined in Section [II Bl 



A. Stress exerted on the clamps 

Long polymer chains have a Gaussian distribution, be- 
coming anisotropic if they contain mesogenic molecules. 
The elastic free energy density of a nematic rubber in 
response to a deformation gradient tensor A can be ex- 
pressed as [l| 

1 / t i \ 1 / Det ( l ) \ 

where fi = n s k^,T is the shear modulus in the isotropic 
state (n s is the number of network strands per unit vol- 
ume). The tensors £ and £ Q generalize the Flory step 
length, whence directions parallel and perpendicular to 
the nematic director n have values £\\ and £±. The ma- 
trix I describes the current Gaussian distribution after 
deformation A, while £ Q gives the initial step lengths. As 
rubber changes shape at constant volume Det (A) = 1. 
Taking n along the z-axis, I assumes the diagonal form 
I = Diag (l±,£x>£\\)' We adopt a simple freely jointed 
polymer model, with a step length a in the isotropic state. 
The elements of I are t\\ = a(l + 2Q) and £± — a(l — Q), 
with Q being the nematic order parameter. Although 
crude, this model quite accurately describes a wide range 
of main-chain and side-chain LCEs [l|, [l3[ . We are only 
concerned with derivatives of F with respect to A and we 
shall suppress A-independent terms in Eq. ((1|). 

Consider an elastomer in an initial nematic state at 
temperature T, with | fl = Diag (i° ± , . Illumina- 

tion changes £ Q to £^ = Diag ($±,£±,£^j- A free sam- 
ple suffers uniaxial spontaneous photodeformation A^ = 

Diag (l/-i/Ap, 1/ WXp, A p ) directed along n = z, with its 
principal contraction A p < 1, and perpendicular elonga- 
tion 1/ •^/A p " due to incompressibility Det (A p ) = 1- The 
free energy density ([TJ corresponding to A = A p is 

Minimisation over A p gives the spontaneous photocon- 
traction A ps = (£\£^ /^|j) 1/3 . Clearly, A ps < 1 
since £fJ£j_ < £^/£°i, that is the sample becomes less 
anisotropic on illumination. 

Clamping (see Fig. [1]), in effect stretches the sam- 
ple by A = 1/A ps along n to restore the sample to the 
length before illumination. It is known [l3j that strain 
little perturbs the underlying nematic order. We as- 
sume that £ p is unchanged after this stretching. Taking 

A = A-A = Diag (1/a/AA ps , l/VAA ps , AA ps ) in © and 
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FIG. 1. A schematic of a nematic elastomer with regular 
mesogenic molecules (white rods) and photo-actives (gray). 
In the dark all photo-active molecules are found in their lin- 
ear trans state (left panel). Upon illumination some photo- 
active molecules undergo a transition to the bent cis state 
(right panel), which would lead to contraction if the elastomer 
were freely suspended. Such contraction is prevented by the 
clamps, resulting in a force / on the clamps. 



the above value for A ps , one gets 




(3) 



The result is the same as that of a classical elastomer 
with a renormalized shear modulus p. If the formation 
state were isotropic followed by cooling to a nematic be- 
fore illumination, F is as in §3§ with a different non-light 
dependent prefactor [![. 

We denote the area of the sample in the nematic state 
before illumination by Aq and its length along the di- 
rector by Lq. Incompressibility gives AqLq = A ps L ps , 
where A ps and L ps are the area and the length of a free 
sample after illumination. The force exerted by a photo- 
elastomer due to the stretching A is 



, _ , . dF _ dF 

J — AqLq — — Apg-j^;-, 



d(XL P 



dX 



(4) 



where A ps = Aq/X ps . After differentiation, we take 
A = l/Ap S since there is clamping. The stress a = f /Aq 
exerted on the clamps takes the form 



a = p 




(5) 



As the shape of the sample is not changing upon the 
illumination, true and engineering stresses are the same. 



The step lengths before illumination are £jj = a(l + 
2Q(T)) and i\ = a(l - Q(T)). Upon illumination, the 
cis isomers act as impurities which lower the nematic 
- isotropic transition temperature, which can also be 
seen as a light-dependent increase of the temperature 
T 0. The elements of £ are B = o(l + 2Q(T cfi )) and 

C\ = a(l— Q(T c g)), where T e g is a fictitious effective tem- 
perature which is the actual temperature T increased by 
a light-dependent term AT, whence T cff = T + AT. We 
estimate order change under illumination by taking shift 
along the Q(T) relation of the dark state. We can exam- 
ine the effect of cis impurities on the free energy. The 
nematic mean field potential is U(9) = —JQP2(cos(8)), 
where 9 is the angle a rod-like mesogen makes with the 
nematic director, J is the nematic mean field coupling 
constant, and Pi is the second-order Legendre polyno- 
mial. The coupling constant J depends on the concen- 
tration p of linear rods as p 2 [ij. We assume for sim- 
plicity that bent cis isomers have no nematic order; they 
weaken the effective potential experienced by the linear 
rods. One can show that the temperature-dependent part 
of Landau - de Gennes expansion of the Maier-Saupe [lH 
free energy has the form Fl = \Aq{T—T*)Q 2 + . . ., where 
T* = J/(5/cb) is the transition temperature. Upon illu- 
mination, one can view the change of p, and thus change 
of J and T* , as an effective change in T at constant T*. 

Our nematic photo-elastomers contain regular meso- 
genic molecules and azo dyes that contribute to the ne- 
matic order when in the trans state. Let po denote the 
total concentration of all mesogenic molecules in the dark 
state (no cis molecules present), and S the molar fraction 
of azo dyes. Upon illumination, the total concentration 
of linear molecules after time t is p(t) = ph + pt(t), where 
Ph = Po(l — S) is the concentration of regular mesogenic 
molecules (constant in time), and pt is that of azo dyes 
in the trans state at time t (pt(t = 0) = poS). Dye 
molecules in trans and cis states contribute to the to- 
tal dye molecules concentration: poS = pt(t) + p c (t), 
where p c (t) is the concentration of cis molecules at time 
t. Expressed in terms of trans n t = pt/po$ and cis 
n c = Pc/pod number fractions, the previous relation be- 
comes nt{t)+n c (t) — 1. For the total concentration of lin- 
ear molecules p{t) one can write p(t) = po[l — S+n t (t)5] = 
po[l — n c (t)S]. Now, the coupling constant J becomes 
J(t) = J [l - n c {t)6} 2 « J [l - 2n c {t)5] for n c (t)6 < 1. 
Thus, cis impurities decrease T* by AT = 2T*n c S, 
which is equivalent to increasing the real temperature 
to T off = T + AT = T + 2T*n c S. 

The actual order parameter (upon illumination) at 
temperature T is then the dark state order parame- 
ter Q(T) shifted to Q(T off ) = Q(T + AT) = Q(T + 
2T*n c 6). If T is not too close to the transition tem- 
perature, T n j, the last expression can be approximated 
by [6( Q(T e g) Q{T) + bn c , where we have introduced 
b = 2(dQ/dT)T*S < 0, since (dQ/dT) < 0. Lineariza- 
tion leads to a simple stress off-dynamics a n (t) = e~*/ r 
(to be discussed in Section lll B[) which is compatible with 
experimental findings. The exact expression for Q(T c ff) 
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(without linearization) converts the simple exponential 
time decay of n c to a non-exponential time response of 
a n . On substituting A ps in Eq. (|5|) and writing the el- 
ements of step length tensors £ and £ in terms of the 

=0 =p 

order parameters Q{T) and Q(T c g) = Q(T) + bn c , one 
gets a cumbersome expression for a. Keeping only lin- 
ear terms in bn Cl consistently with linear decay processes, 
one obtains 



-36/i 



(1-Q(T))(1 + 2Q(T))' 



(6) 



The stress a is of course positive since b is negative. The 
cis number fraction is time and depth dependent n c — 
n c (x,t), which leads to a (x,t) dependence of the stress. 
To calculate the average stress a(t), that is a measure 
of the force, exerted on the clamps one must integrate 
over depth a(t) — a{x, t)dx/d, where d is the sample 
thickness. To do this, we have to explore the behavior of 
the cis number fraction n c (x,t). 



B. Light absorption 

The light intensity I varies with depth x (see Fig. Q} 
due to absorption by trans and cis species of dye 
molecules 

dl 

— = -jT t I(x,t)nt(x,t) - jT c I(x,t)n c {x,t), (7) 

where T t and T c are rate coefficients for photon absorp- 
tion. Here 7 = HujpqS subsumes the energy fuv of each 
absorption of a photon from the beam and the absolute 
number density of chromophores pqS. To simplify the 
above relation we normalize I(x,t) by the incident inten- 
sity Iq to obtain an X = I/Iq. The combinations 7l\ and 
7r c will be written as l/rf t and l/d c respectively. Now 
Eq. ((7]) becomes 



dX 

dx 



n t (x,t) 



l(x,t) 



(x,t) 



T{x,i). 



(8) 



Assuming that the the irons-population of absorbers 
does not change appreciably, n t (x,t) ~ 1, one obtains 
Beer attenuation X = e~ x / dt . To close Eq. (jH]) one needs 
the rate equation for the irons-population at (x,t), 

= -r) t T t I(x,t)nt(x,t) + (i] c T c I(x 7 t) + - J n c (x,t). 

The changes in n t are due to photoconversions (with 
quantum efficiencies nt per photon absorption of trans 
to cis transition, and rj c per photon absorption of cis to 
trans transition), and thermal back reaction from the cis- 
population at a rate 1/t, with r being the cis lifetime. 
Eq. ^ can be rewritten as 



dn t 
' dt 



(1 + /3X(x, t))-[l + (a + /3)X(x, t)]n t (x, t) , (10) 



where the combinations a = ri t r t IoT and /3 = j] c T c Iot 
measure how intense the incident beam is compared 
with the material constants It = l/faT^r) and I c = 
1/(i] c T c t). The parameter a is the ratio between the for- 
ward trans to cis conversion rate 77tr t /o) and the thermal 
back rate, 1/t; parameter /3 has similar interpretation. 
Note that both a and (3 depend on temperature T since 
the cis to trans thermal decay is activated, r = t(T), 
and t is measurable directly or in stress decay. 
After integration of Eq. ((8]) over depth x we get 



n c {x, t)dx 



d - d t A{d, t) 
1 - d t /d c 



(11) 



where the absorption is determined by the relation 
A(x, t) = — In 1(2;, t). From the above equation one eas- 
ily obtains the average stress 



a(t) 



3b[i 



\-{d t /d)A{d,t) 



(1-Q)(1 + 2Q) 1- dt/dc 



(12) 



In experiments one usually measures the normalized 
stress a n (t) — W(t)/a s , where W s is the stationary state 
stress, 



cfn(t) 



l-{d t /d)A(d,t) 
l-(d t /d)A{d) ■ 



(13) 



Here we adopt the convention that A(d) denotes the 
steady state value of the absorption A(d, t) . The same 
convention will be used for X, nt and n c . 

For off-dynamics, after setting J = in Eq. (J9|), 
one gets a simple equation for the cis number fraction 
n c (x,t) = n c (x, 0)e _t / T , where n c (x,0) represents the 
spatial profile of n c at the instant the light is switched off 
t = 0. On inserting n c (x,t) into Eq. (j6|), averaging over 
depth, and normalizing the average stress with its value 
at t = 0, one gets a simple expression a n (t) — e~ t / T . 

Given that the stress W n (t) for on-dynamics depends on 
both A(d, t) and A(d), we shall firstly discuss the behav- 
ior of these quantities. In the steady state, dn t /dt = 0, 
we have 



n t (x) 



1 + /3X{x) 
1 + (a + f3)X(x) 



, n c (x) 



aX{x) 



1 + (a + p)l{x) 

(14) 

Now the equilibrium solution of Eq. (J8]) can be expressed 
in the form 



ln(I) 



a — T]/3 

W+v) 



In 



l + /3(l + n)X 
1 + P(1 + V) 



; (15) 



where n = r/t/r) c is the ratio of the quantum efficien- 
cies. The above expression for X{x) provides a general- 
ization of the usual Beer's law. Deviations from Beer's 
law come about because at high intensities the cis popu- 
lation increases and is generally less absorbing than the 
trans species. In the limit T c — > 0, the parameter /3 — > 
so that Eq. (fTS"]) reduces to 



ln(Z) + a(X - 1) 



x 
"d+? 



(16) 
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FIG. 2. Light intensity versus the reduced depth x/dt for 
various reduced times t/r. The dash-dotted lines represent 
the case with the presence of back photoconversion (a = 36, 
P — 1, and r\ — 3 in this example), while the solid lines 
correspond to the absence of back photoconversion (a = 36, 
P = 0); note that even a small value of P (compared to a) 
significantly affects the variation of I with depth. In the limit 
of large t/r one achieves the stationary regime described by 
Eq. (T5J. 



For a — 0, one finds the standard Beer's law X = e 
while for large a the law acquires the linear form 



(17) 



at least over depths up to x ~ adt whereupon X is small 
and the ln(2T) again prevails to give a finally exponen- 
tial penetration. The variation of light intensity with 
reduced depth for large a is represented in Fig. [5J by the 
solid line t/r — oo in the case of absence of back photo- 
conversion {(i — 0), and by the dash-dotted line t/r = oo 
when /3 ^ 0. Note that even a moderate value of (3 (see 
Fig. [2]) has a major impact on the variation of X with 
depth. Non-Beer absorption was first explored for dyes 
in nematic liquid crystals by Statman and Janossy [9j 
and by Corbett and Warner [liJ-Qjj], and experimentally 
by Serra and Terentjev (l5j . 

Beer absorption has no dynamics since it holds only 
if the number fraction n t is unchanging, n t — 1. To 
explore the dynamics of non-Beer absorption, we have to 
solve the coupled Eqs. © and (jTUJ) . Indeed, using these 
two equations one gets |12| 



OA 
~dt 



-A+ — + (a + /3)(X-l) + — + 



I 
dt 



Xdx. 

Clearly, the absorption A can be expressed in terms of 
reduced variables x = x/d t and t — t/r. Taking deriva- 
tive of the above equation with respect to x, and using 



FIG. 3. The cis number fraction n c against the reduced depth 
x/dt for various reduced times t/r. The dash-dotted lines 
represent the case a = 36, P = 1, and r\ — 3 while the solid 
lines represent the case a = 36, P = 0. In the large time limit 
one approaches the solution presented by the second equation 

of col. 



X — e , we get a partial differential equation in A(x, t), 

^= = 1 - ^ [1 + (« + I3)e- A ] + 13(1 + r,)e~ A . (19) 

We analyze this equation numerically using the boundary 
condition A(0,t) = and the initial condition A(x, 0) = 
x. The results are presented in Fig. [2] As expected the 
light intensity decreases with reduced depth more rapidly 
in the presence of back photoconversion (f3 =t 0, dash- 
dotted lines) than in the absence of back photoconver- 
sion (f3 — 0, solid lines). Note that in both cases the 
stationary regime is reached for moderate values of the 
reduced time t/r. 

Before proceeding with the analysis of the behavior 
of the stress, it is tempting to consider briefly the time 
evolution of the cis number fraction n c = n c (x, t). Thus, 
by taking I from Eq. (fTU)) and substituting it into Eq. © , 
one gets 



(n G + n c ) 



1 - Hr. 1 



1 - Hr. 1 



VP 



= [h' c (n c - 1) - <(1 + ho)] + - (h' c n c - h c n' c ) , (20) 

a 

where for simplicity we used the notation h c = dn c /dt 
and n' c = dn c /dx. In the numerical integration of this 
equation we use the boundary condition 



nc(0j) 



l + a + P 



{l-exp[-t(l + a + P)]}, (21) 



and the initial condition n c (x, 0) = 0. The quoted bound- 
ary condition is obtained by integration of Eq. (|10p for 



6 




FIG. 4. Approximate thickness Xc/dt of cis layer against the 
reduced time t/r for a — 36, f3 = (solid line) and a = 36, 
P = 1, and = 3 (dash-dotted). Both curves reach their 
stationary values after a moderate time (measured in units of 
r) which coincide with corresponding values (thin solid lines) 
obtained from Eq. (|22p as explained in the main text. 



x = 0. Figure [3] shows function of the reduced 

depth x/dt for various reduced times t/r. It seems that 
the stationary solution, given by the second equation 
of (THf . is reached already for moderate values of the 
reduced time t/r. 

We estimate the thickness x c /dt of cis layer at t/r by 
calculating positions of points of the curves from Fig. [3] 
where second derivatives of n c over x/dt change the sign. 
In the stationary case, this condition can be expressed 
analytically, by taking second derivative of the second 
equation of dl~4l) with respect to x/dt and using (|8]), 

(a + 0)0(1 + 7y)I c 2 + 2(a - rj/3)l c -1 = 0, (22) 

where I c = I(x c /dt). Inserting the positive solution of 
this equation into Eq. (fl"5j) we obtain an equation for 
the stationary value x c /d t . In the non-stationary case 
we proceed numerically. Results for a = 36, /3 — 0, and 
a = 36, /3 = 1 are presented in Fig. 2) In both cases x c /dt 
increases approximately logarithmically with time in the 
interval from about 0.5r to 1.5r; this can be easily seen 
if one presents curves from Fig. |4]on a logarithmic scale. 
The stationary solution is reached in about t = 4t for 
(3 = 1 and t = 5r for /3 = 0. It is also interesting to com- 
pare the time needed to reach the stationary state with 
finite /? with the time needed to convert equally thick 
layer without cis absorption. In Fig.S]this ratio is about 
10, so even though the cis absorption is small, it strongly 
affects the photodynamics of the system. The back pho- 
toconversion has also been found to be important in the 
analysis of experimental findings [l6j|. 



III. RESULTS 

As we have seen, the absorption A at the back of the 
sample, x = d, in Eq. (TT3]) for stress is actually a func- 
tion of reduced variables, A(d, t). Then one can generate 
cf n (i) for different values of parameters a, /3, d and r. Ex- 
perimental data for a n (t) was fitted Q to the simple em- 
pirical form a n (t) — 1 — exp[— (t/r on ) ,3on ], with /3 on < 1. 
This stretched exponential form must fail at short times, 
and fitting at long times is difficult. We have fitted our 
theoretical a n (t) to the above form. In the absence of 
cis to trans photoconversion (/3 = 0) we can get fits for 
/Son > 1 only. If we allow, however, back photoconversion 
we can obtain agreement with the stretched exponential 
form, /3 on < 1. It is, therefore, very important to take 
into account back photoconversion to successfully fit the 
experimental data. 

We consider the experimental data for compound 
SCEAzo2-c-10 in Fig. 5(a) of Ref. 0. There are several 
quantities entering our relation (|13|) for stress: the ratio 
of the quantum efficiencies 77, the characteristic thermal 
relaxation time r, thickness of the sample d — d/dt in 
units of the characteristic length, and the dimensionless 
parameters a and /3. For the ratio of the quantum effi- 
ciencies we adopt the estimate 77 « 3 of Ref. [17]; some 
larger values of 77 (up to 77 « 4) have also been reported. 
Our analysis shows, however, that the quality of the fit 
is not very much affected by the particular value of 77 in 
the range 77 w 3 — 4. For the thermal relaxation times 
t(T) we take the experimental values found in the stress 
off-dynamics 0: r = 51.5, 33, 23, 15, 10.4, 7.2 min, for 
temperatures T = 45, 50, 55, 60, 65, 70 °C respectively. 
With this choice the number of required fit parameters 
is reduced to three, d, a, and /3. 

We first fit the experimental data for a n (t) at To = 
45 °C to our expression (fl~3"|) . using d, a and /3 as fit 
parameters; see Fig. [5] The optimal values of fit param- 
eters were found to be a = 54.9, (3 = 1.7, and d = 20. 
Given that the sample thickness [8| was d = 300 /j,m, the 
corresponding Beer length is dt — 15 /im, which is not 
far from an independent estimate dt ~ 10 71m obtained 
from the azobenzene absorption spectra. This rough es- 
timate of dt can be obtained by assuming that the mo- 
lar extinction for SCEAzo2-c-10 at its absorption maxi- 
mum coincides with the corresponding known value for 
the azobenzene in benzene [181 ]. Let us note, however, 
that quite good fits can be also obtained for some other 
values of fit parameters. This ambiguity can be settled 
obviously by a further reduction of the number of the fit 
parameters - for example by measuring d. Anyway, in 
our case situation is not so serious, because the parame- 
ter d is temperature independent, while parameters a and 
/3 depend on temperature T only through thermal relax- 
ation time, a(T) = 77 t r t / r(T) and /3(T) = 77 c r c / r(T). 
For example, at other temperatures T, one can write 
a(T) = a(T )T(T)/ T (To) and /3(T) = (3(T )t(T)/t(T ). 
Thus taking the experimental values of thermal relax- 
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FIG. 5. Normalized stress as a function of time. The- 
oretical predictions are presented by solid lines, while ex- 
perimental data, taken from f|j, are presented by different 
symbols corresponding to different temperatures. For bet- 
ter visibility the five curves corresponding to temperatures 
T > 45 °C are shifted along the a n -axis by the values 
0.05, 0.1 , 0.15, 0.2, 0.25 respectively. Only the T = 45 °C 
curve is fitted, while all other curves are theoretical. 



IV. CONCLUSIONS 

Wc have demonstrated that both significant force mag- 
nitude and complex force dynamics result from nonlin- 
ear optical absorption. A bleaching wave of increased cis 
concentration, and hence contribution to retractile force, 
penetrates a sample with a highly characteristic dynam- 
ics in which the (small) absorption of the cis moiety is 
essential. Fitting W n (t) at one temperature yields the rel- 
evant material parameters of the dye responsible for the 
opto-mechanical response. The a n (t) response at other 
temperatures is then obtained with no further fit by sim- 
ply scaling these parameters by the separately measured 
decay rates at the other temperatures. This astonishing 
predictive power points to the validity of the nonlinear 
temporal-spatial optical absorption model used. 

Future experiments should separately measure Beer 
penetration depths in the weak absorption limit, and 
the nonlinear material parameters. In that event there 
should be no fit parameters at all, as in this current work 
away from the reference temperature. With the material 
parameters thus measured, one could employ the theory 
to examine more complex systems such as non-clamped 
elastomers and systems where the response is more com- 
plex due to director patterning. 



ation times t(T) for temperatures quoted in Fig. [3] we 
determine the corresponding values of a(T) and P(T). 
Then using these values and expression (|13|) for the nor- 
malized stress we simply plot corresponding curves for 
temperatures T > Tq without any fitting. Note excellent 
agreement of these theoretical predictions with experi- 
mental data. Extracting a and /3 by fitting at one To 
gives universal, fit-free agreement at other temperatures. 
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